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Abstract 

Entanglement within a given device provides a potential resource for quantum information processing. Entanglement between 
system and environment leads to decoherence (thus suppressing non-classical features within the system) but also opens up 
a route to robust and universal control. The latter is related to thermodynamic equilibrium, a generic behavior of bi-partite 
quantum systems. Fingerprints of this equilibrium behavior (including relaxation and stability) show up already far from the 
thermodynamic limit, where a complete solution of the underlying Schrodinger dynamics of the total system is still feasible. 
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1. Introduction 

There have been various attempts to reduce thermo- 
dynamics to some underlying more fundamental the- 
ory. While the vast majority of the pertinent work done 
in this field has been based on classical mechanics [1] , 
a reduction to quantum mechanics has also attracted 
increasing interest [2] . Decoherence [3] has, during the 
last years, often been discussed as one of the main 
obstacles for the implementation of large-scale quan- 
tum computation. Quantum thermodynamics [4], on 
the other hand, tries to show that decoherence is far 
from being just a technical nuisance but a generic phe- 
nomenon of partite quantum systems giving rise to 
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some of the most dominating, if classical, features of 
closed (finite) quantum systems: thermal equilibrium. 

The program of a quantum foundation of thermody- 
namics should contain the following points: 
(1) a definition of thermodynamic quantities (in the 
ideal case, as a function of microstates); (2) a deriva- 
tion of the second law of thermodynamics under appro- 
priate constraints (including stability, irreversibility, 
universality); (3) a justification of the Gibbsian funda- 
mental form (state functions and conjugate variables); 
(4) a proof of extensivity or intensivity, respectively, 
of the thermodynamic variables; (5) a characterization 
of thermodynamic systems (as opposed to other sys- 
tems) ; (6) a kind of correspondence principle (explain- 
ing the efficiency of standard classical approaches de- 
spite their underlying quantum nature). 

Here we cannot do justice to this rather challenging 
program. Instead we want to address some of the main 
results available to us. 
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2. The Model 



3.1. Accessible Region (AR) 



We consider a bipartite system - an observed system 
or gas g and an environment or container c, described 
by the Hamiltonian 

H = H S + H C + 1 . (1) 

Note that also the environment requires a full quantum 
treatment and should not be replaced by (classical) 
boundary conditions. I defines the interaction between 
these two subsystems. 

Weak coupling between system and environment has 
routinely to be assumed in standard thermodynam- 
ics [5]; otherwise the concept of intensive and exten- 
sive variables would lose its meaning. Furthermore, in 
this case the full spectrum of the coupled system will 
not look significantly different from the one that re- 
sults from a mere convolution of the two spectra of the 
uncoupled system. To quantify the weak coupling pre- 
condition, we require 

\f(P) « (H s ), (H c ). (2) 

This inequality must hold for all states that the to- 
tal system can possibly evolve into under given con- 
straints. 

The weak coupling has further to be classified. Not 
so much for practical, but for theoretical reasons, the 
most important contact conditions are the microcanon- 
ical and the canonical conditions. In the microcanoni- 
cal contact scenario no energy transfer between system 
and environment is allowed, as opposed to the canoni- 
cal contact. 



3. Microcanonical Conditions 

If a system is thermally isolated, it is not necessarily 
isolated in the microscopic sense, i.e., not uncoupled to 
any other system. The only constraint is that the in- 
teraction with the environment should not give rise to 
any energy exchange. As will be seen later, this does 
not mean that such an interaction has no effect on the 
considered system, a fact that might seem counterin- 
tuitive from a classical point of view. This constraint, 
however, leads to an immense reduction of the region 
in Hilbert space which the wave vector is allowed to 
enter. This reduced area is called "accessible region" 
of the system. 



If the energies contained in the gas g and the envi- 
ronment c, respectively, 

E s := (H s ) , E c := (H c ) (3) 

are to be conserved, i.e. if these two energies are con- 
stants of motion, the following commutator relations 
should hold 

[HZ,H)=0, [H c ,H}=0. (4) 
It then follows from 

[H s , H] = [H s , H s ] + [H s , H c ] + [H s , /] = (5) 
that 

[hM]=q, [H c J}=0. (6) 

Except for these constraints we need not specify I in 
more detail. All interactions that fulfill this relation will 
create perfectly microcanonical situations, regardless 
of their strength or any other feature. And, as will be 
shown, there are a lot of possible interactions that do 
fulfill these conditions and create entanglement and 
therefore give rise to the increase of local entropy. 

Due to (4) the local energy projectors P\ of the gas 
system and P B of the container 

PZ = 52\A,a){A,a\, P B = £|B, b)(B, b\ (7) 

a b 

commute with the full Hamiltonian, 

[PZ,H] = [H,H]=0. (8) 

Here a(b) specify the degenerate states (degeneracies 
Na(Nb)) for given energy eigenvalue A(B) in the gas 
(container) system. Thus, because the system is not 
allowed to exchange energy with the environment the 
joint probability W AB must be conserved 

wpxPbm = |v>rw| 2 = e l^wf 

a,b a,b 

= W AB , (9) 

and is set by the initial state. This means that the en- 
ergy probability distribution {W AB } is a constant of 
motion. Vice versa, any state that features this same 
energy probability distribution as the initial state be- 
longs to the accessible region and could possibly be 
reached during microcanonical dynamics. 
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In the following we mainly consider initial product 
states, states that have zero local entropy in the begin- 
ning and for which 

£ |vd s (o)| 2 = k« (o)f k s (o)| 2 = w A w B . 



a,b 



a,b 



(10) 

This is the only constraint that microcanonical condi- 
tions impose on the accessible region of Hilbert space. 
Note that this does not mean that local entropy is con- 
stant. 

3.2. The "Landscape" of P 9 in the Accessible Region 

To demonstrate that the largest part of the accessible 
region is filled with states of almost minimum purity 
(maximum entropy), we proceed as follows: 

(i) First we compute the (unique) state with 
the lowest possible purity, p^ in (with purity 
P(Pmin) ~ ^"min) triat i s consistent with the 
given initial state and the microcanonical con- 
ditions, consequently with a given energy prob- 
ability distribution {W^}. 

(ii) Then we compute the average of P g over the 
total accessible Hilbert space region, denoted by 
lP e l 

(iii) We show that this average purity is very close 
to the purity of the lowest possible purity state 
Pmin f° r a large class of systems. Considering 
only these systems, which then define the class of 
thermodynamic systems, we can conclude that 
pg pj, pj\ n for almost all states within the ac- 
cessible region. Note that this conclusion is only 
possible because of the fact that the purity of 
Pmin i s the absolute minimal purity which can 
be reached at all in the system. A quantity with 
a mean value close to a boundary cannot vary 
very much. Thus it is not possible that the dis- 
tribution of the purity within the accessible re- 
gion is something else but a very flat "lowland" , 
with a "soft depression" at p^ in (see Fig. 1) and 
a "peak" with P g = 1. 

(iv) Since all states from the accessible region have 
the same energy probability distribution {W A } 
(remember (10)) and the minimum purity state 
Pmin i s consistent with this distribution, all 
other states within the accessible region that 
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Fig. 1. Qualitative picture of the purity landscape in the 
microcanonical case. The biggest part of the accessible re- 
gion is at P P^ in or at P = ■ There is however only 
a small zone featuring P significantly above P^ in or at the 
extreme P = 1. The only topological property this rough 
picture refers to is the relative size of different regions. 

states that are very close to p^ in (in this context 
close means in terms of the distance measure 
Tr {(p 9 - pLJ 2 }) - Thus, as long as the trajec- 
tory keeps wandering through the compartment 
filled with those states the gas system is locally 
in a stationary state, equilibrium is reached. 
Details can be found in Ref. [4,6]. 

3.3. Microcanonical Equilibrium 

The minimal purity for subsystem g under the given 
constraints is 



P s . =y 

min ^ 



(11) 



For the Hilbert space average of P s in the accessible 
region we find 

W. 



v N * v u 

B \ A 

r2 u/2 I 



+ E 



A,B 



WXW 2 B {N A +N B ) 
N A N B + 1 



(12) 



feature P s 



must yield reduced local 



This average is thus a unique function of the invariants 
W A , W B , specified by the initial product state and the 
degeneracies N A , N B . 

If the degeneracy of the occupied energy levels is 
large enough so that 
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N A N B + 1 » N A N B , (13) 

which should hold true for typical thermodynamic sys- 
tems, (12) reduces to 

A A B s 

The first sum in this expression is obviously exactly 
-Pmin' so that f° r systems and initial conditions, in 
which the second sum is very small, the allowed region 
almost entirely consists of states for which P g « P^ in . 
The second sum will be small if the container system 
occupies highly degenerate states typical for thermo- 
dynamic systems, in which the surrounding is much 
larger than the considered system. This is the set of 
cases mentioned already in Sect. 3.2: all systems fulfill- 
ing this pre-condition are called now thermodynamic 
systems. Thus we can conclude that all states within 
the accessible region are very close to p^ in and have 
approximately the purity P^ in . The density operator, 
which has P s = P s . and S s = SB,,., , and which is con- 

mm iiidJt } 

sistent with the microcanonical conditions, is unique. 
The density operators with P s m P^ in should not devi- 
ate much from this one and should therefore also have 
S s « Sm ax , the latter being 




4. Energy Exchange Conditions 

4.1. The Accessible and the Dominant Region 

Our approach to the "energy exchange conditions" 
will be based on similar techniques as before. The pos- 
sibility of a partition is still assumed. But now there is 
no further constraint on the interaction I, since energy 
is allowed to flow from one subsystem to the other. The 
only constraint for the accessible region therefore de- 
rives from the initial state of the full system, and the 
fact that the probability to find the total system at 
some energy E, 

W{E):= Y, W *b=Y, ( 16 ) 

A,B/E A,B/E a,b 

should be conserved, where A, B/E stands for: all A, B 
such that E A + E° B = E. This constraint is nothing 
else but the overall energy conservation. 




Dominant Region 



Fig. 2. Qualitative picture of the purity landscape. In the 
canonical case the accessible region contains a dominant re- 
gion which almost entirely fills the accessible region. Within 
the dominant region, all states feature the same energy 
probability distribution. Thus all topology from the micro- 
canonical case (cf. Fig. 1) transfers to the dominant region. 

One could try to repeat the above calculation under 
the energy conservation constraint, but now it turns 
out that the average purity over the accessible region 
is no longer close to the actual minimum purity. Fur- 
thermore, the energy probability distribution of the in- 
dividual considered system is no longer a constant of 
motion. Thus, we proceed in a slightly different way: 

(i) Contrary to the microcanonical case the prob- 
ability to find the gas (container) subsystem at 
some given energy is no longer a constant of mo- 
tion here. But one can prove that there is a pre- 
dominant distribution, {M^g}, which almost all 
states within the allowed region have in common. 
The subregion formed by these states is called 
the "dominant region" . 

(ii) Having identified the "dominant region" we 
demonstrate that this region is by far the biggest 
subregion in the accessible region of the system 
(see Fig. 2). 

(iii) Once the existence of such a dominant region has 
been established, we can use the results from the 
microcanonical conditions to argue that almost 
all states within this dominant region, which is 
specified by a fixed energy probability distribu- 
tion for the considered system, feature the max- 
imum local entropy that is consistent with the 
predominant distribution. Based on this analy- 
sis we get the equilibrium state of the considered 
system (see Sect. 4.2). 

Just like in the previous case our subjective lack of 
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knowledge about where to find the system within the 
accessible region is irrelevant. The reduced local state 
p s (t) as a function of the full state \ip{t)) should always 
evolve into a state with a fixed probability distribution 
Wa, and an almost time invariant entropy, which is 
the maximum entropy consistent with this (canonical) 
distribution. Nevertheless, the state of the full system 
continues to move in Hilbert space with the constant 
velocity v (v » energy uncertainty). 

Again, details can be found in Ref. [4,6]. 

4.2. The Equilibrium State 

Finally, to find the marginal, dominant energy prob- 
ability distribution W% of the gas system individually, 
one has to sum the compound probabilities W%b over 
the irrelevant container system to obtain 

N A N B 



wi = *£wi B = Y: 

B B/E 



B/E 



N C {E B )W{E) 
N(E) 



(17) 



where again the sum over B/E denotes a summation 
under the condition E — E A + E B , and N A = N g (E A ) 
and N B = N C (E B ) are the respective degeneracies. 
Since E B = E - E\ is a function of E for fixed E A we 
switch from a summation over B to a summation over 
E, 



(18) 



This is the energy probability distribution for the gas 
system that one will find with overwhelming probabil- 
ity for a thermodynamic system. Simply by exchanging 
the indices (up to here everything is symmetric with 
respect to an exchange of the subsystems) we find the 
marginal dominant energy probability distribution for 
the container system. 

So far we have only established the energy probabil- 
ity distributions for almost all states from the acces- 
sible region, but nothing has been said about entropy, 
purity, etc. The equilibrium state is still undetermined. 
Once the trajectory has entered the dominant region, 
we can assume that the trajectory will practically never 
leave it, because this region fills almost the whole ac- 
cessible region of the system. But since all states within 
the dominant region feature the same energy proba- 
bility distribution, motion within the dominant region 



will never give rise to any further energy exchange be- 
tween the subsystems. As a consequence the situation 
is as if the system was controlled by microcanonical 
conditions. 

Therefore, we can take the arguments from Sect. 3.3 
to identify the equilibrium state. Following this idea, 
we can start with (14) and use the dominant energy 
distribution W AB , finding for the Hilbert space average 
of the purity of the gas 



[^i*E^-+E 



(19) 

A "A B " 

Once more it is possible to conclude that the second 
term (due to the environment) is much smaller than 
the first one for a sufficiently degenerate environment. 
The first term is exactly the minimum purity of the 
gas system within the dominant region. Thus, almost 
all states from the dominant region will yield approx- 
imately the same local gas state too. This equilibrium 
state p c s q is, again, the state of minimum purity (max- 
imum entropy) that is consistent with the dominant 
energy distribution, 



(20) 



One problem remains: the dominant energy proba- 
bility distribution W A (18) is not independent of the 
initial state since different energy probability distri- 
butions of the local initial state may result in different 
overall energy probability distributions W(E), and 
those clearly enter (18) and thus even (20). Normally 
the canonical contact of standard thermodynamics 
leads to an equilibrium state, which does not depend 
on the initial state. The canonical contact turns out 
to be a special subclass, as we will demonstrate in the 
next Section. 



5. Canonical Conditions 

For a canonical situation the gas system should relax 
into the canonical equilibrium state, independent of 
the initial conditions. This behavior can be found, if a 
further condition is taken into account: a special form 
of the degeneracy of the environment N B . 

Let us assume an exponential increase of the con- 
tainer degeneracy 

N B = NSe aE °B, (21) 
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where a, N§ are some constants. Such a degeneracy 
structure is typical for modular systems [4]. We start 
again from (18) using (21) for the degeneracy of the 
environment 

~ aE W(E) 



N 



a e 



N(E) 



(22) 



Obviously, the sum does not depend on A at all. Since 
W% has been constructed as some probability distribu- 
tion it is still normalized by definition. Therefore the 
sum has to reduce to a normalizing factor. Finally we 
get for the dominant energy probability distribution of 
the gas system 



wi = 



N A 



(23) 



This result is no longer dependent on the initial state! 

The energy probability distributions of almost all 
states from the accessible region consistent with the 
constraints is then the canonical distribution: Since the 
argumentation for the minimal purity state (state of 
maximal entropy) remains unchanged, the equilibrium 
state reads now 

1 NT- 



Z A N A e-" E A 



J2e- aE A\A,a){A,a\. (24) 

A, a 



Obviously, this is the well known canonical equilibrium 
state with the inverse temperature f3 — a. 



value for two systems in energy exchanging contact, 
and if the energy of a system is changed without 
changing its volume, it should be a measure for the 
energy change per entropy change. 

Most definitions rely on the second property. 
Maxwell connected the mean kinetic energy of a clas- 
sical particle with temperature. In the canonical en- 
semble (Boltzmann distribution) it is guaranteed that 
the energy change per entropy change equals temper- 
ature. And the ensemble mean of the kinetic energy of 
a particle equals fcsT in this case. Thus, if ergodicity 
is assumed, i.e., if the time average equals the ensem- 
ble average, temperature may indeed be defined as 
the time averaged kinetic energy. Similar approaches 
have been proposed on the basis of the microcanonical 
ensemble [8,9]. However, if temperature was given by 
a time average over an observable, the proper averag- 
ing time remains open and thus the question on what 
minimum timescale temperature may be defined. Fur- 
thermore, this definition was entirely based on ergod- 
icity. Nevertheless, it allows, at least to some extent, 
for an investigation of processes, in which temperature 
varies in time and/or space, since that definition is not 
necessarily restricted to full equilibrium. 

To avoid those problems of standard temperature 
definitions, we present here yet another, entirely quan- 
tum mechanical definition. 



6. Temperature 



7. Definition of Spectral Temperature 



One could claim that temperature should only be 
defined for equilibrium and thus there was no need to 
define it as a function of the micro state. Based on this 
reasoning temperature would then simply be defined 

as 

1 dS d , 1 dG(E) 



k B T dE dE V ' G(E) dE 



(25) 



with G(E) being the state density. In this way one 
would neglect all dynamical aspects (see [7]), since this 
definition is based on the Hamiltonian of the system 
rather than on its state. Strictly speaking, this defini- 
tion would exclude all situations in which temperature 
appeares as a function of time or space, because those 
are non-equilibrium situations. 

A quantity like temperature is essentially deter- 
mined by two properties: It should take on the same 



We define the inverse spectral temperature as [4] 
W + W M S 



1 

M 

E 



w t + w. 



HNi/Ni. 



(26) 

where W t is the probability to find the quantum sys- 
tem at the energy E iy M is the number of the highest 
energy level E M , while the lowest one is labeled E . 
This formula is motivated by the following idea: For a 
two level system it seems plausible to define tempera- 
ture just from the energy probability distribution and 
the degrees of degeneracy as 



W N, 



exp 



E l 



k B T 



(27) 
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The definition (26) results if one groups the energy 
levels of a multi-level system into neighboring pairs, 
to each of which a "temperature" is assigned via the 
above formula, weighted by the average probability for 
each pair to be occupied. This definition obviously de- 
pends only on the energy probability distribution and 
the spectrum of a system. It thus cannot change in time 
for an isolated system, and it is always defined, inde- 
pendent of whether or not the system is in an equilib- 
rium state. Thus there should be many systems or sit- 
uations with such a temperature, which do not exhibit 
thermodynamic properties at all. The latter will only 
show up in equilibrium situations or close to those. 

If the spectrum of a system was very dense and if 
it was possible to describe the energy probability dis- 
tribution, {Wj}, as well as the degrees of degeneracy, 
{Ni}, by smooth continuous functions (W(E),N(E)) 
with a well defined derivative, (26) could be approxi- 
mated by 
1 



/ 

Jo 



(28) 

This can further be simplified by integrating the first 
term to yield 

^-L «W(0)-W(£ max ) 



+ 



W(E) dN(E) 



dE . 



(29) 



N(E) dE 

Since for larger systems typically neither the lowest nor 
the highest energy level is occupied with considerable 
probability (if the spectra are finite at all), it is the 
last term on the right hand side of (29) , that basically 
matters. This term can be interpreted as the average 
over the standard, system based, rather than micro 
state based definition of the inverse temperature. 



AE 
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Fig. 3. Microcanonical scenario: A non-degenerate 
two-level-system (gas) is weakly coupled to a system with 
one energy level of degeneracy N c = 50. This is a model for 
a system in contact with a much larger environment such 
that no energy can be exchanged between the system and 
environment. 

condition of energy conservation. What sort of physi- 
cal structure could give rise to those spectra is not con- 
sidered here, since it turns out to be irrelevant to some 
extent. 

The spectra are translated into discrete diagonal ma- 
trix Hamiltonians that describe the decoupled bipar- 
tite system. These can be chosen to be diagonal, with- 
out any loss of generality, since any such system may 
be analyzed in the energy eigenbasis of its decoupled 
parts. The form of the interaction depends on the con- 
crete physical subsystems and their interactions. But 
since the "guess" is that for the quantities considered 
here (entropy, occupation probabilities, etc.) the con- 
crete form of the interaction should not matter, the in- 
teraction is taken as some random matrix thus avoid- 
ing any bias. The interaction has to be "small" . In this 
way it is hoped that we get models that are as "typ- 
ical" for general thermodynamic situations as possi- 
ble. Many of the situations analyzed in this Chapter 
are very similar to those treated within the context of 
quantum master equations. But note that in order to 
apply the theories at hand to those systems, neither a 
Markovian nor a Born assumption has to hold. We sim- 
ply solve the respective Schrodinger equation in finite 
dimensional space. 

8.1. Entropy under Microcanonical Conditions 



8. Equilibrium Properties of Model Systems 

We now turn to some numerical data, based on a 
certain type of models. These models are still rather 
abstract and may thus be viewed as models for a whole 
class of bipartite systems. The subsystems g, c are spec- 
ified by their respective spectra only, or, rather, those 
respective parts which play any role at all under the 



All data in our first example refer to a situation de- 
picted in Fig. 3 (cf. [10]). The "gas" (the system under 
consideration) consists of a two-level system, both lev- 
els being non-degenerate (N§ = Nf — 1), while the 
"container" (the environment) consists of just one en- 
ergy level with degeneracy A?" c = 50. This is necessarily 
a microcanonical situation regardless of the interaction 
I. The container cannot absorb any energy, therefore 
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Fig. 4. Relative size of Hilbcrt space compartments: this 
histogram shows the relative frequency of states with a 
given local entropy S, among all states from the accessi- 
ble region. In this case the maximum possible entropy is 
■Smax = 0.423 /cb ■ Obviously, almost all states feature en- 
tropies close to the maximum. 

energy cannot be exchanged between the systems. In 
this situation the probabilities to find the gas system 
in the ground (excited) state W§ (Wf) are conserved 
quantities and in this example chosen as 

W§ = 0.15 , Wf = 0.85 . (30) 

As described in Sect. 3, the Hilbert space average of the 
purity of the gas system is under condiction (13) given 
according to (14) by JP g J = 0.765. The corresponding 
minimum purity (9) is P„ lin = 0.745. As explained in 
Sect. 3, we find here 

I^ g ] « PL, , (31) 

a situation, in which almost the entire accessible region 
would be filled with the compartment containing only 
states of almost maximum local entropy. 

To confirm this expectation, a set of random states, 
uniformly distributed over the accessible region, has 
been generated. Their local entropies have been calcu- 
lated and sorted into a histogram. Since those states 
are distributed uniformly over the accessible region, the 
number of states in any "entropy bin" reflects the rela- 
tive size of the respective Hilbert space compartment. 

The histogram is shown in Fig. 4. The maximum lo- 
cal entropy in this case is = 0.423 fca- Obviously, 
almost all states have local entropies close to S^^. 
Thus compartments corresponding to entropies of, say, 
S s > 0.4 fee indeed fill almost the entire accessible re- 
gion, just as theory predicts. Local pure states (S e = 
0) are practically of measure zero. 



0.45 




50 100 150 200 250 300 

Fig. 5. Evolution of the local entropy for different initial 
states. A universal state of maximum entropy (equilibrium) 
is reached, independent of the initial state. 



In order to examine the dynamics, a coupling I is 
introduced. To keep the concrete example as general 
as possible, I has been chosen as a random matrix in 
the basis of the energy eigenstates of the uncoupled 
system, with Gaussian distributed real and imaginary 
parts of the matrix elements of zero mean and a stan- 
dard deviation of 

AI = 0.01A.E . (32) 

This coupling is weak, compared to the Hamiltonian 
of the uncoupled system. Therefore the respective in- 
teraction cannot contain much energy. The spectrum 
of the system (see Fig. 3) does not change significantly 
due to the coupling, and afterall the environment is not 
able to absorbe energy. 

Now the Schrodinger equation for this system, in- 
cluding a realization of the interaction, has been solved 
for initial states consistent with (30). Then the local 
entropy at each time has been calculated, thus get- 
ting a picture of the entropy evolution. The result is 
shown in Fig. 5. Obviously the entropy approaches 
S'max within a reasonable time, regardless of the con- 
crete initial state. Thus the tendency towards equilib- 
rium is obvious. The concrete form of the interaction 
/ only influences the details of this evolution, the equi- 
librium value is always the same. If the interaction is 
chosen to be weaker, the time scale on which equilib- 
rium is reached gets longer, but, eventually the same 
maximum entropy will be reached in any case. 
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Fig. 6. Canonical scenario: A two-level gas system is weakly 
coupled to a three level environment, such that energy can 
be exchanged. The exponential degeneracy scheme of the 
container system guarantees a full independence of the equi- 
librium state from the initial state. 

8.2. Occupation Probabilities under Canonical 
Conditions 



□00 



Fig. 7. Evolution of the ground level occupation probability 
for three different random interactions. The dotted line 
corresponds to a weaker interaction. Even in this case the 
same equilibrium value, = |, is approached, only on a 
longer timescale. 



The second model analyzed numerically is depicted 
in Fig. 6. The considered (gas) system, again, consists 
only of a non-degenerate two-level system. The envi- 
ronment (container) in this case is a three-level system 
with an exponential "state density" : N B — 50 • 2 s with 
B — 0, 1, 2. This has been chosen since theory predicts 
for such a degeneracy scheme of the environment an 
equilibrium state of the gas system, which should be 
independent of its initial state (see (22)). If we restrict 
ourselves to initial states featuring arbitrary states for 
the gas system but container states that only occupy 
the intermediate energy level, no other container levels 
except for those given could be reached, even if they 
were present: This is due to energy conservation and 
holds for the limit of weak interactions /. 

In this case the model can also be seen to represent a 
situation with a much larger environment and we find 
from (23) 

W d {El)= 2 -, W A {Ef) = \. (33) 

To keep the situation as general as possible, / was, like 
in Sect. 8.1, chosen to be a matrix with random Gaus- 
sian distributed entries in the basis of the eigenstates 
of the uncoupled system, but now with energy transfer 
allowed between the subsystems. 

For this system the Schrodinger equation has been 
solved and the evolution of the probability to find the 
gas system in its ground state, W(Eq) is displayed 
in Fig. 7. The different curves correspond to different 
interaction strengths, given by the standard deviation 



of the distribution of the matrix elements of I, AI: 

A/ so ii d . daBh ed = 0.0075AS , A/ dottcd = 0.002AS . 

(34) 

Obviously, the equilibrium value of W d (i?Q) = 2/3 is 
reached independently of the concrete interaction I. 
Within the weak coupling limit the interaction strength 
only influences the timescale on which equilibrium is 
reached. 

Figure 8 displays the evolution of the same probabil- 
ity, W (Eq) , but now for different initial states, featur- 
ing different probabilities for the groundstate, as can 
be seen in the figure at t = 0. The equilibrium value is 
reached for any such evolution, regardless of the spe- 
cial initial state, thus we confirm the effective attractor 
behavior typical for thermodynamics. 

Figure 9 shows the evolution of the local entropy of 
the gas system for the same three initial states as used 
for Fig. 8. 

The maximum entropy, consistent with the equi- 
librium value of the energy probabilities, is S^ ax = 
0.637 ks- This is also the value one finds, if one max- 
imizes entropy for fixed mean energy (Jaynes' princi- 
ple [11]). Obviously, this value is reached for any ini- 
tial state during the concrete dynamics of this model. 
This supports the validity of (20), which states that 
the density matrix of the equilibrium state is diagonal 
in the basis of the local energy eigenstates. 

To analyze the formation of a full Boltzmann distri- 
bution, we turn to the system depicted in Fig. 10. Here 
the "gas" system is a non-degenerate equidistant five- 
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Fig. 8. Evolution of the ground level occupation probability 
for different initial states. The theoretically predicted equi- 
librium value is reached, independent of the initial states, 
as expected for canonical conditions. 



Fig. 9. Evolution of the local entropy for different initial 
states. S = 0.637 &b is the maximum entropy that is consis- 
tent with the equilibrium energy probabilities. This maxi- 
mum entropy state is reached in all cases. 

level system the container system a five-level system 
with degeneracies N% = 6 ■ 2 s (B = 0, . . . , 4), which 
should lead to a Boltzmann distribution. We restrict 
ourselves to initial states, where for both subsystems 
only the intermediate energy level is occupied (symbol- 
ized by the black dots in Fig. 10). Due to energy con- 
servation other states of the container system would 
not play any role in this case even if they were present, 
just like in the previous model. Figure 11 shows the 
probabilities of the different energy levels to be occu- 
pied W(E\). While the gas system starts in the inter- 
mediate (third) energy level, soon a Boltzmann distri- 
bution develops: Obviously, each probability becomes 
twice as high as the one for the level above. This is ex- 
actly what theory predicts (see (18)) for the environ- 
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Fig. 10. Canonical multi-level scenario: A five-level gas sys- 
tem is weakly coupled to a five level container system with 
an exponential degeneracy scheme, such that energy may 
be exchanged. Black dots symbolize the initial state. This 
set-up should lead to a Boltzmann distribution. 




Fig. 11. Evolution of the energy occupation probabilities. 
After some relaxation time a Boltzmann distribution is 
reached. Each probability is twice as high as the one for 
the next higher energy level, as theory predicts. 

ment degeneracy scheme under consideration. 



8.3. Probability Fluctuations 

Finally, we consider fluctuations in time due to finite 
size effects imposed by the environment. For this pur- 
pose a system almost like the one depicted in Fig. 6 is 
analyzed, but now with a degeneracy scheme given by 

N% = ^ • 2 s . (35) 

The ratios between the degrees of degeneracy of the 
different container levels are thus the same as for the 
system sketched in Fig. 6, but the overall size of the 
container system is tunable by Nf. For various Nf, 
the Schrodinger equation has been solved numerically, 
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Fig. 12. Fluctuations of the probability for the considered 
system to be in the ground state AjWj, in dependence 
on the number of eigenstates of the environment system 
iVf . Obviously, the fluctuations go down, with increasing 
environment. 

and the following measure of the fluctuations of the 
occupation probability of the ground level of the gas 
system has been computed 



A 2 t W§ := 



1 



tf - U 



for initial states with 

W E (0) = 0.2 



(W ( f(t)fdt 



W§(t)dt 



WftO) = 0.8 



(36) 



(37) 



Figure 12 shows the dependence of the size of these 
fluctuations on the container system size iVf . The small 
crosses are the computed data points, the dashed line 
is a least square fit to a function proportional to 1/iVf , 



AtWi = 



0.053 

JVjf 



(38) 



This fit is in very good agreement with the theoretical 
expectation, although the trajectories are not ergodic. 



9. Summary and Conclusions 

The reduction of thermodynamics to an underlying 
theory has challenged physicists for more than a cen- 
tury. So far such attempts have met with partial suc- 
cess only. In this context an interesting and, as we be- 
lieve, consistent route is offered by quantum thermo- 
dynamics. 



Quantum thermodynamics as a concept has not yet 
a generally accepted meaning (see, e.g. [12] for a dif- 
ferent opinion). Our version has been introduced as a 
special field within decoherence theory. It makes use 
of a special methodology, which is based on statisti- 
cal techniques to investigate the structure of a prod- 
uct Hilbert space under certain constraints. These con- 
straints, though, are by themselves emergent proper- 
ties of the system environment model. 

Within this quantum thermodynamics we are able 
to show that thermodynamics is already contained in 
the deterministic Schrodinger dynamics. Fingerprints 
of equilibrium behavior show up already in suprisingly 
small bipartite quantum systems. As a consequence 
thermodynamics need not be considered as a final 
remedy after all detailed microscopic descriptions 
have failed: Deterministic Schrodinger dynamics of 
the whole and (apparent) relaxation behavior of the 
embedded part is no contradiction at all. 

Such a picture should also give new insight on ther- 
modynamic applications in the nano-regime [13]. This 
expectation is corroborated by the fact that pertinent 
models then have to be on the quantum level anyway, 
while the conditions for a thermodynamic description 
can be made explicit, as this is part of the complete 
quantum description of system and environment [14] . 
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